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ABSTRACT 

-  - i  A  systeai  of  two  alliptie  p.d.e.'s,  which  Model  the  conduction-convection 
problem  in  a  porous  Medium  with  change  of  phase  is  considered. 

The  first  equation  describes  the  heat  conduction  and  holds  in  a  fixed 
domain.  The  second  takes  into  account  the  convective  Motions  and  holds  in  the 
unknown  Melted  part  of  the  region. 

The  existence  of  a  locally  regular  weak  solution  is  proved  by  Making  use 
of  various  compactness  arguments. 
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SIGNIFICANCE  AND  EXPLANATION 


Consider  the  problem  of  artificial  freezing  or  thawing  of  ground  water. 
The  problem  has  important  engineering  applications  such  as  ground 
stabilization  at  construction  sites  and  the  prevention  of  water  leakage  into 
tunnels  or  shafts. 

In  the  melted  part  of  the  ground  the  convective  motions  forced  by 
temperature  gradients  affect  significantly  the  free  boundary  of  separation 
between  the  melted  and  solid  regions. 

As  a  mathematical  model  of  this  situation  tie  study  a  system  of  two 
elliptic  equations,  one  describing  heat  conduction  and  the  other  the  velocity 
flow  of  convection.  We  demonstrate  that  the  problem  has  a  solution  in  a 
suitable  weak  sense.  Regularity  properties  of  the  solution  are  also  derived. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


EXISTENCE  FOR  A  PROBLEM  IN  GROUND  FREEZING 


E.  Di  Benedetto* ' 1  and  C.  M.  Elliott2 

1.  INTRODUCTION 

There  axiata  laboratory  and  thaoratical  avidanca  to  suggest  that  convaction  can  play 
an  iaportant  rola  in  haat  transfer  with  change  of  phase  in  a  porous  medium  (see 
16,7,8,1]).  We  mention  specifically  the  artificial  freezing  or  thawing  of  groundwater  in 
view  of  the  important  engineering  applications  such  as  ground  stabilisation  at  construction 
sites  and  the  prevention  of  water  leakage  into  tunnels  or  shafts  (see  [11,  12]).  In 
certain  circumstances,  the  moving  boundary  problem  considered  by  the  above  authors  has  a 
nontrivial  steady  state.  This  leads  to  the  study  of  a  free  boundary  problem  for  a  system 
of  two  elliptic  partial  differential  equations. 

Let  fl  be  a  bounded  domain  in  R2  lying  between  two  closed  curves  3+fl  and  3_Q. 

The  domain  SI  is  the  union  of  two  sets  SI  and  S)*  which  are  respectively  the  frozen  and 
unfrozen  regions  with  the  phase  change  temperature  being  zero  at  the  unknown  free  boundary 
T  separating  the  sets.  Figure  1  can  be  interpreted  as  a  vertical  cross-section  of  a 
portion  S)  of  frozen  ground,  where  a  horizontal  pipe  of  cross-section  3+S)  carries  warm 
oil.  Then  around  the  pipe  there  is  an  unfrozen  region  Sl+.  If  (2+  and  S)~  are 
Interchanged  then  one  could  regard  Figure  1  as  depicting  the  artificial  freezing  of  a 
horizontal  layer  of  ground  by  the  insertion  of  a  vertical  freeze  pipe.  In  this  case  the 
experiments  of  Frivlk  and  CoaU.nl  [6]  and  the  ntsserical  work  of  Goldstein  and  Reid  [7]  and 
Barrett  and  Elliott  [1],  show  that  the  flow  of  groundwater  can  significantly  affect  the 
extent  of  the  frozen  region 


^Department  of  Mathematics,  Indiana  Uhiversity,  Bloomington,  IN  47405 
2 Department  of  Mathematics,  Imperial  College  of  Science  and  Technology,  Huxley  Building, 
Queen's  Oats,  London  SW7  2BZ 


*  Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041.  Partially 
supported  by  National  Science  Foundation  Grant  48-296-80  MCS  8300293. 


Si 


Figure  1 

Let  u, v  and  p  denote,  respectively,  dimensionless  temperature,  velocity  and 

pressure.  The  model  equations  are  then  the  steady  state  convection-conduction  equation, 

the  continuity  equation  and  Darcy's  law  which  links  the  pressure  to  the  velocity  t1,7). 

The  problem  can  be  formulated  as  follows. 

♦ 

Find  {u.v.p.r}  such  that 

*  2 

(1.1)  -Au  ♦  Xv*Vu  -  f.  +  £  t —  f  in  SI 

°  i«1  “xi  1 

d.2)  s)4  =  tx  e  $)|u(x)  >  oil  si"  s  {x  e  si|u(x)  <  o) 

a  -  si4  u  r  u  a~»  ssi  s  a+a  u  j_a 

(1.3)  v  ■  -kQ{Vp  +■  g(u,x)}  in  S)+ 

(1.4)  div  v  «■  0  in  S)+;  v  =  0  in  0 

(1.5) 

(1.6)  u  “  ♦  on  3S1,  ♦  >  0  on  3+8,  ♦  <  0  on  3_0  . 

2  3 

Here  X  >  0  is  a  prescribed  physical  constant,  f.  ♦  l  f.  represents  a  known 

0  i-1  3xi  1 

distribution  of  heat  sources  (or  sinks)  in  SI,  kQ  is  the  coefficient  of  permeability  and 
+ 

q  is  a  buoyancy  force. 

Remarks:  (i)  If  f ^ ,  i  *  1,2  is  the  Heaviside  function  in  the  x^  direction  then  on 
the  right  hand  side  of  (1.1)  we  allow  sources  (or  sinks),  like  Dirac  masses. 
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£  -  , 
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(li)  The  problem  has  been  foraulated  having  in  Bind  the  case  when  thawing  occurs  and 
0  is  a  vertical  layer.  This  is  then  a  problea  in  free  convection  where  the  flow  is 
driven  by  the  vertical  temperature  gradient.  In  the  case  of  artificial  freezing  in  a 
horizontal  layer  the  effect  of  a  prescribed  flow  on  the  extent  of  the  frosen  region  is  of 
interest  and  then  (1.S)  nay  be  replaced  by  the  prescription  of  a  non-zero  normal  velocity 
or  pressure. 

We  observe  that  the  incompressibility  condition  (1.4)  and  the  Darcy  law  (1.3)  imply 
that  the  pressure  solves  the  elliptic  equation 

(1.7)  div  kQ(Vp  ♦  g(  u,x) ]  •  0  in  Q+  . 

The  aim  of  this  paper  is  to  prove  an  existence  theorem  for  (1.1)-(1.6)  in  a  suitable 
weak  sense.  The  difficulty  here  is  represented  by  the  fact  that  (1.3)  and  (1.7)  are 
conditions  which  hold  on  the  unknown  region  St* .  They  might  be  interpreted  in  the  sense  of 
distributions  over  fl*  if  8*  is  open.  This  tatter  information  will  be  implied  by  the 
local  HSlder  continuity  of  the  temperature.  The  Min  points  of  the  proof  are  the 
introduction  of  a  penalized  problem  (in  the  spirit  of  (3,4])  and  a  careful  use  of  the 
intonation  (1.4)  and  the  fact  that  the  space  dimension  is  two,  in  order  to  prove  that  u 
is  HSlder  continuous.  We  now  come  to  a  precise  formulation  of  our  results. 

Definition.  By  a  weak  solution  of  (1.1)-(1.7)  we  mean  a  triple  (u,v,p)  such  that 
(i)  u  e  w1,2(Q)  n  i,“(Q)  n  ca(0)»  u  -  ♦  on  80 
(11)  v  e  (L2(Q)]2  and  v  «  0  a. a.  on  (u  <  0) t  v*H.  _  “  0 
(Hi)  p  e  W,,2(tu  >  on, 
satisfying  (the  simulation  notation  is  used) 

(1.8)  /  (Vu  -  XvujVndx  -  /  f.ndx  -  /  f  n  dx 

a  a  0  a  1  *i 

wn  e  w1,2(0)  n  l“(Q) 

(1.  9)  /  [Vp  ♦  g(u,x)]«Vndx  “  0 

lu>0] 

vn  e  w1,2(0)  such  that  sup  n  c  (u  >  0)  , 


-3- 


(1.10) 


/  (v  -  k.[Vp  ♦  g(u,x)])*ndx  -  0 
[u>0] 

Vn  e  [L2(9)]2  such  that  supp  n  C  [u  >  0]  . 
a 

Remarks!  (1)  Since  we  require  u  e  C  (0),  the  set  (u  >  0)  is  open  in  the  relative 
topology  of  9  and  hence  (1.9)-(1.10)  are  meaningful. 

(li)  (1.9)— (1.10)  imply  that  div  v  -  0  in  V ( (u  >  0]). 

The  integrals  in  (1.8)-(1.t0)  are  trail  defined,  modulo  basic  assiaptions  listed  below. 
[Aq]  39  satisfies  the  cone  condition  (see  [10]) 

2+K 

[A,]  f±  e  L  (9),  K  >  0,  i  -  0,1,2 
1  /2  • 

(AjJ  ♦  e  H  '  (39)  n  l  (39),  ♦  >  0  on  3+9  and  +  <  0  on  3_9 

♦  3  2 

[Aj]  g  i  R  ♦  R  is  continuous  and 

|g(u,x)|  <  Clul*1  +  gr0< >t)  , 

where  C  and  t  are  given  nonnegative  constants  with  t  e  [0,1),  and 

gQ(x)  e  Lj(9)  . 

Theorem.  Imt  [a01-[A3]  hold.  Then  problem  (1.1)-(1.7)  has  a  weak  solution. 

Remark.  let  t  »  1  in  [Aj) .  The  existence  of  a  weak  solution  in  this  case  follows  from  a 

remark  in  section  5  provided  the  data  is  sufficiently  small . 

The  paper  is  organized  as  follows.  In  section  2  a  family  of  penalized  (or  smoothed) 

problems  is  introduced.  Basic  results  for  this  family  are  stated  in  section  2  and  proved 

in  sections  4  and  5.  The  theorem  is  demonstrated  in  section  3.  Standard  notation  for 

function  spaces  is  employed  throughout  the  paper  and  1*1  is  used  to  denote  the  LP(Q) 

P«9 

norm.  The  measure  of  a  set  A  in  iP  is  denoted  by  |A|. 


2.  CONSTRUCTION  OP  THE  SOLUTION 


For  e  >  0  consider  the  Lipechitz  continuous  function  s  ♦  ke(s)  defined  by 


(2.1) 


k£(s) 


k0  *  6  (k0  "  6,1 


8  >  0 

s  e  [-e,0] 

c  s  <  -e 

♦ 

end  the  problem  of  finding  functions  u£ ,v£ ,p  satisfying 

u£  e  tf*'2(Q)i  ue  "  ^  on  ^ 


(2.2) 

(2.3) 

(2.4) 

and 

(2.5) 


v£  8  (L2(0)]2»  div  v£  •  0  in  D'(0) 


pe  €  W1,2(Q) 


/  (Vu  *Vn  -  Xw  u  Vn}dx  -  /  f  ndx  -  /  f.n  dx,  vn  e  w1,2<n) n  l“(0) 
ne  Ee  n°  a  1  xi 

/  k«- *  **.- *  lvP*  +  9(uc#*)J *Vndx  •  0,  »n  e  w1,2{a> 
fl  e  e  e  e 

/  (v£  -  k£lue>fvPe  +  g(ue,x)})ndx  *  o,  vn  e  [l2ui)}2  . 


nils  smoothing  of  the  problem  (1.1)-(1.2)  has  a  sound  physical  basis.  In  reality  it 
may  Indeed  be  the  case  that  the  coefficient  of  permeability  takes  a  tiny  non-zero  value  in 
the  frozen  region  (see  (61).  Hence  the  study  of  (2.1)— (2.5)  is  of  interest  in  its  own 
right. 

The  proof  of  the  theorem  rests  upon  the  following  propositions 
Proposition  1. 

For  all  e  >  0,  problem  (2.2)— (2.5)  has  a  solution.  Moreover,  Ve  >  0 


(2.6) 

(2.7) 

(2.8) 


'V-.O  *  G0 


'V2,n  ♦  ,Vue'2,0  4  G1 


,/Vu>  I’Pe'S.a  «  G2 


-5 


(2.9) 


upon  the  data  and  Independent  of  e  such  that 

(2. to)  |uc(x)  -  U£(y)|  <  r(K)|x  -  yl“.  V(x.y)  0  K,  PC  >  0  . 

Remark 

If  the  boundary  datua  ♦  e  C°,  then  the  equi-H51der  continuity  of  the  net  (u£) 
carries  up  to  Q,  with  constants  y  and  a  depending  upon  the  date  and  a. 


The  proof  of  propositions  1  end  2  are  postponed  to  sections  4  and  5. 
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3.  PROOF  or  TOE  THEOREM 

Let  us  assume  Propositions  1,2  for  tha  moment  and  let  ua  conclude  the  proof  of 
existence  by  uaing  thaaa  facta. 

Because  of  (2.6)-(2.10)  aubnata  out  of  {v£>  and  (u£)  can  ba  selected  (and 
relabeled  with  e)  such  that 

v£  ♦  v  weakly  in  lL2(Q)l2 
u£  ♦  u  weakly  In 

u£  ♦  u  uniformly  on  compact  aeta  of  (2  . 

Letting  e  ♦  0  in  (2.3)  proves  (1.8).  From  (2.4)-(2.5)  we  have 

/  v  •  Vndx  -  o  me  w1,2(a) 
a 

and  hence 

div  v  “  0  in  D' (0)  . 

Since  (u£)  are  equibounded  and  equi-HSlder  continuous  in  Q,  the  uniform  limit  u 
is  H5lder  continuous  in  fl,  and  therefore  the  sate  [u  >  0)  and  (u  <  0)  are  open  in  the 
relative  topology  of  fl. 

Let  K  be  compact  and  contained  in  (u  <  0) .  there  exists  o  >  0  and  eQ  <  o/4 
such  that 

u(x)  <  -o  ,  ft  6  It 

u£(x)  <  -c/2  ,  vx  e  It  ,  Ve  <  eQ 

ke(ue>  -  c  ,  ft  6  it  ,  ve  <  eQ  . 

By  (2.5) 

v£  -  -ke(u£)(7pe  ♦  g(u£,x)]  a.e.  in  Q 

and  therefore 

/  |v  |2dx  <2  /  k2(u  )|Vp  |2dx  *  2  f  k2(u  )|g|2dx 
K  6  K  *  6  6  K  C  6 

<  2  el/kc(ue)  IVpe!«|  Q  ♦  2  e2i;(ue,*).2fQ 
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Letting  c  ♦  0  by  the  weak  lower  seaicontlnuity  of  the  norm  we  obtain  v  »  0  a.e.  In 

K,  and  aince  K  la  an  arbitrary  compact  set  in  [u  <  0]  we  have 

♦ 

v  «  0  a.e.  In  tu  <  0]  . 

It  remains  to  prove  ( 1 .9)-( 1 . 10) .  Let  K  c  [u  >  0)  be  compact.  There  exist 
o  >  0,  Eg  >  0  such  that 

u(x)  >  a  ,  Vx  e  K 

ue(x)  >  o/2  ,  Vx  6  K  ,  VC  <  eq  . 

k£(uc(x))  “  kp  ,  Vx  e  K  ,  VE  <  E0  . 

By  (2.8)  and  the  definition  of  k£(u£) 

'7Pe'2,K  «  C2//r0  * 

The  vectors  {Vp£}  are  therefore  weakly  compact  in  [L2 (K ) ] 2  and  for  a  subnet,  relabeled 
with  ek 

Vpr  *  K  weakly  in  [l2(K)12  . 

EK  K 

Let  (K  }  be  a  family  of  compact  subsets  of  [u  >  0]  such  that 
n 

K  C  [u  >  0],  K  C  K  §  u  K  -  (U  >  0]  e 
n  n  nt-i  .  n 

n>i 

Then  by  a  diagonal  process  a  subnet  out  of  {Vp£}  can  be  selected  (and  relabeled  with 
E )  such  that 

7PC  ♦  t  weakly  in  (L2(  [u  >  0])I2 

Next  we  want  to  Identify  t  as  the  gradient  of  a  function  p  e  w''2([u  >  0)). 

N  _  «• 

Let  A  be  an  open  set  in  R  ,  N  >  2  and  let  Po(A)  be  the  space  of  all  CQ  vector 
fields  ♦  in  A  which  are  solenoidal  (div  $  •  0).  He  set 

J  (A)  -  (the  (L2(A)]N  closure  of  P  (A)}  . 
a  o 

Let  also  6(A)  denote  the  space  of  those  square  integrable  vector  fields  in 
[L2(A)]**  obtained  as  gradients  of  functions  in  W1,2(a). 
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By  Weyl's  1. 


mu,  [9],  the  orthogonal  decomposition, 

[L2(A)JN  -  J  (A)  •  G(A)  , 
o 

holds . 

He  must  show  that  {  e  G([u  >  0]),  or  equivalently  that 

(3.1)  /  !•$  dx  -  0  ,  e  J  (tu  >  0]) 

[u>0]  0 

Since  (u  >0])  is  dense  in  J^([u  >0])  it  will  suffice  to  prove  (3#1)  for  all 

l  e  Vgll U  >  on. 

How  Ve  >  0  e  w1,2((u  >01),  therefore  for  all  |  e  8  ([u  >  0])  we  have 
€  0 

/  Vp  dx  «  /  p  div  {  dx  «  0  . 

Iu>0]  (u>0 J 

Letting  e  ♦  0  (3.1)  follows.  He  will  set 

t  -  Vp,  p  e  w1,2([u  >  on  . 

Consider  now  (2.4)  and  (2.5)  where  n  and  T|  are  supported  on  same  compact  set 

K  C  [u  >  0].  By  an  argument  similar  to  the  one  above,  if  e  is  sufficiently  small  we  must 

have  ke(uc)  -  k0.  Letting  e  ♦  0  we  recover  (1.9)  Vn  e  w1,2(n)  supp  n  C  K,  and  (1.10) 
■»  2  2 

vn  e  [L  (Q)l  supp  n  C  K.  Since  K  C  (u  >  0)  is  arbitrary  (1.9)  and  (1.10)  follow. 


4.  A  PRIORI  ESTIMATES 

It  is  the  purpose  of  this  section  to  prove  the  following  two  la— as. 

"Hrrs  1 

Let  v  e  (L2(8)]2,  div  v  -  0  in  P'(8)  and  Ivl  <  V  for  a  given  positive 

•l" 

constant  V.  If  (Agl-tAjl  hold  then  the  problem 


(4.1) 


f  J  (Vu  -  Avu}Vn  dx  -  /  ffln  dx  -  /  f  n  dx  vn  e  w1,2(8)  l~{8) 

a  q  a  *i 


has  a  unique  solution  u  e  W1,2(8).  Also  the  following  a  priori  estimates  hold: 

(4.2,  'u'j.n  *  <  V1  +  v2) 

(4.3)  ,u,-.a  <  V1  +  v) 

and  u  in  HSlder  continuous  on  compact  subsets  of  8  with  the  HBlder  exponent  a 
and  Holder  constant  Y  depending  on  [AqJ-IAj]  and  V. 


m 

Let  u  e  L  (8)  be  arbitrary  and  C  a  given  positive  constant.  If  [AQ;A31  hold 

1  2 

then  the  problem  of  finding  p  “  p(u,c)  e  W  *  (8)  such  that 

(4.4)  /  k  (u)[Vp  ♦  g(u,x)]-Vndx  -  0  Vn  e  W1,2(8) 

a 

1  2 

has  a  solution  irtiich  belong  to  a  unique  class  [p]  e  W  '  (8)/K.  The  vector  fie'd 
v(u,£)  “  k£ (u) [Vp(u,£ )  ♦  g(u,x>]  is  unique  and  the  following  estima tes  hold: 

(4.5)  l/ie(uj  |Vp(u,e)|l2>Q  <  k’/2{c|8|1/2lul^fl  ♦  »90«2#n) 

(4.6)  i;<a*e)l2.S  <  21t0{C,fl,,/2,u,-,8  *  'V2.81  ’ 


where  C,  t  and  g«  are  introduced  in  [A.].  Finally  if  u  ,u.  e  L  (8)  then 


(4.7)  IVptu^C)  -  7P<U2'e,l2,B  *  G1l,U1  "  U2,»,Q  *  ,l«*tt1»x>  "  9<U2'X>  * ‘-.Q1 


*  «'V2.B  *  'Vi.0  *  'V-V 


lv(ut,e)  -  »<“2*e,,2,Q  *  G2,,P(u|'e)  “  Vp<u2'C>l2,Q 


where  depend  on  |Q|,  k„,  C  end  e.  In  perticuler  as  e  ♦  0  the  beccae 


unbounded. 


4-(a)  Existence,  uniqueness  and  boundedness  for  u 


The  available  existence  and  uniqueness  theory  cannot  be  directly  applied  here  since 


v  Is  not  regular  enooqh.  Let  us  prove  uniqueness  first.  Suppose  u,  and  u2  are 
solutions  of  (4.1).  Then  w  ■  u,  -  Uj  *  %# 1 ' 2 (0)  satisfies 


/  (T«  -  Xvw)7ndx  -  o  »n  e  w1,2(0)  n  l"(Q)  . 
Q 


Taking  n  -  w  «  w1,2(Q)  In  (4.»),  where 

H 


n  ,  w  (x)  >  n 
fl 


Wn*X>  “  I  e  l*n#nl 


-n  ,  w(x)  <  -n 


and  noting  that  dlv  v  »  0,  results  In  the  equation 


/  |Vw  |2dx  -  0  n  -  1,2,... 


Thus  the  convergence  of  w  to  w  inplies  the  equation 


/  | Vw| 2dx  -  0 

a 


and  the  uniqueness  of  u. 


Existence  is  proved  by  considering  a  sequence  of  probleas 


(4.10)  u  -  ♦  e  W,,2(0)  :  /  (Vu  -  Xv  u  )Vndx  -  /  ( f nn  -  f.n  )dx  Vn  e  W1,2(0) 
n  q  n  n  n  q  u  1 


i 


where  4  la  an  extension  of  4  into  3  such  that 


141  ,  ,  <  4„ 

m1,2(B)  0 


for  a  given  fixed  constant  4^,  and  (vn>  is  a  sequence  in  1)^(3)  which  converges  to 


in  J  (8)  such  that 
0 


•v  I,  „  <  Ivl,  _  =  V  . 
n  2,3  2,3 


Since  div  v  »  0,  Poincare's  inequality  implies  that 
n 


( (w,n) )  5  /  (7w  -  Xv  w)Vndx 
3  " 


is  an  inner  product  on  M1,2(S)  which  qeneratea  an  equivalent  topology  In  W,,2(B). 
Similarly  the  right  hand  side  of  (4.10)  can  be  written  as  f(n)  where  f(«)  is  a 
continuous  linear  functional  on  M1,2(S)  and,  since  4  C  W1,2(8),  a  standard  application 
of  the  Riesz  representation  theorem  yields  the  existence  and  uniqueness  of  a  solution  to 
(4.10). 

Me  now  obtain  some  estimates  for  which  are  independent  of  n.  Taking 

n  “  u  -  4  in  (4. 10)  we  obtain 


(4.11)  /  (|V(u  -  4)|2  -  ~  Av  *V(u  -  4)2}dx 

_  n  4,  n  n 


/  {Av  V(U  -  4)4  -  V4V(u  -  4)}dx  ♦  /  {f.(u 

4  pi  n  4  n  n 


‘,  +  ,i^  ,Un  -  •»* 


Observing  that 

/  v  *V(u  -  4)2dx  -  0  , 

3  n  n 

2  2 

a  repeated  application  of  the  inequality  ab  «  aa  ♦  b  /o,  yields  the  existence  of  a 
constant  Y  depending  only  on  A  such  that 
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9 


<♦ 


(4.12) 


■*"v  -  •"!.«  <  ’('"-.o'K'U  •  '«'2!.o  j,  ”,,2,a) 


*  ""'Via  *  «■“.  -  "l,:, 


Poincari's  inequality  and  (4.12)  inply  the  estiaute 
(4.13) 


“n'lfl  *  ,7un'2.a  *  V1  *  y2)  * 


where  Yg  depend*  on  the  data  and  Q. 


*n  L  bound  for  u^  will  now  be  derived.  Let  i  >  1$!^  Then  the  function 

f  12 

1  “  (un  -  i)  e  W  *  (0)  and  can  be  used  as  a  test  function  in  (4.10)  to  give 


(4.14)  /  | V(u  -  i)+|2  +  Xv  *7u  (u  -  t)+dx  -  /  f  (u  -  i)  +  -  f  —  (u  -  i)+« 

o  »  nnn  'On  i  ax,  n 


dx 


Setting 


and  noting  that 


=  {«  e  Q|u(x)  >  t) 


/  *  2,+<1*  *  2  /  \‘V<\  -  «> 


+2 


dx 


we  obtain  front  (4.14)  the  inequality 


/  |V(«n-i)+|2dx<.(un-t)+.2#n  +  iIo/  f±2dx 


Holder's  inequality  and  (A^J  imply  that 
(4.15)  -  **+<2  '  -*+*2 


1-2/2+k 


,7(»n-‘)  ,2.«<  '(“„*»>  '2.fl  +  ^l'V 

where  Yj  depends  only  upon  the  data.  Let  M  be  a  number  larger  than  21*1^  3Q  and 
conelder  the  sequence  of  increasing  levels 

-  M  ♦  N(1  -  2”*)  a  -  0,1,2,...  . 

A  bound  on  ess  sup  u(x)  will  be  obtained  by  taking  l  -  l  in  (4.15).  Holder's 
xefl  ■ 

inequality  and  Sobolev  eabedding  (see  for  example  inequality  (2.12)  on  p.  45  of  I10J)  iaply 


13- 


that  for  2  <  p  <  • 


/  -  t^J^dx  <  1/  (un  -  t-)*pdx}2/p|At  |1_2/p 


■+ 1 


<  c|n|1/pIV(u  -  t  >V  q|A  |1_2/p 

n  »+1  2,Q 


where  c  depends  only  on  p.  He  deduce  froa  (4.15)  that 


(4-,6»  '  <“»  -  V/2**  4  72,Av/"2/P,'(U„  -  Vl,+2,2.fl  *  I^J1"272**}  • 


4 

Next  observe  that  since  t. .  «  t  ♦  M/2  , 

m+ 1  i 


-  ^1>+,2,a  4 


and 


-  v*'m  > ;  <■.  -  >  (-*?>’  ■*,  i 

Ag  2  1 

■+1 


Choose 


and  set 


P  -  4(2  +  k)/k  e  (2,-) 


a  -  k/2(2  +  k)  e  (0,1)  . 

Substituting  into  (4.16),  we  obtain,  for  a  constant  depending  only  upon  the  data 


*+1  2c  _irM  2+20 


(4.1?)  /  V(V-)  +(V-)  H/  (»n  - 


Setting 


y-  "  “2  l  <Un  -  *.,+2<Ix 


M  ft 


(4.17)  iaplies 

(4.18) 


v,  < 
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where  Yj  depends  only  upon  the  dets  end  in  perticuler  Is  independent  of  m.  Inequality 
(4.18)  implies,  with  the  aid  of  lemma  4.7  of  (101 ,  p.  66,  that  y  ♦  0  as  *  ♦  “  if 


Y0  <  (V'VV,/ 


Consequently  if  N  is  chosen  to  satisfy 


we  have  that  y_  ■  0  i.e.  /  (u  -  2M)+2dx  »  0  and 


ess^sup  u(x>  <  T5  Q.»nn»2<n)  . 

where  depends  only  upon  the  data  and  is  independent  of  v.  A  bound  for  ess  inf  u(x) 

a 

is  obtained  in  a  similar  fashion  and  we  deduce  that 


(4*,9>  'V-.a  *  *6  "“{,*,..9a',un,2.a)  * 

Passing  to  the  limit  in  (4.13)  and  (4.19)  as  n  tends  to  infinity  we  obtain  the  existence 
of  e  solution  to  (4.1)  setisfying  (4.2)  and  (4.3). 


4-(b)  Hglder  continuity  of  u 

The  results  of  (10]  Imply  that  the  HBlder  exponent  and  constant  of  u  may  be 

estimated  in  terms  of  the  u’(B),  q  >  2,  norm  of  v  and  the  data.  However  the  proof  of 

the  theorem  requires  that  these  estimates  are  independent  of  e  and  it  is  known  only  that 

•vlj  0  <  V  where  V  does  not  depend  on  e,  (see  Section  4-(c>).  In  the  following,  the 

arguments  of  (10,  p.  81,  90]  are  adapted  to  demonstrate  that  the  extra  information  of 
♦ 

div  v  -  0  and  the  apace  dimension  being  two  is  sufficient  to  conclude  that  u  is  HBlder 
continuous  in  Q  with  estimates  depending  only  on  Ivl  and  the  data. 

2pU 

bet  Xq  be  an  arbitrary  point  of  B  and  let  B(p)  denote  the  ball 
B(p)  =  (x  s  |x  -  xQ|  <  p).  For  x0  e  B,  p  will  always  be  chosen  so  small  that 
B(2p)  fl.  For  o  e  (0,1),  c(«)  will  always  denote  a  smooth  cutoff  function  such  that 
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U)  aupp  c  C  B(«) 

(II)  C  =  1  on  B(p  -  op)  and  C(x)  8  10,1} 

(III)  |V?|  <  (op)"1. 

The  Hdldar  continuity  of  u  la  shown  by  tha  judicious  uaa  of  an  Integral  Inequality 
which  follows  directly  from  the  Integral  Identity  (4.1).  For  arbitrary  k  C  I,  the 
results  of  4-(a)  iaply  that  the  functions 


(4.20) 


r\  -  ±(u  -  k)*;2  8  51,2(a> 


are  adniaaahle  teat  functions  in  (4-1)  and  aetting 


p  =  {x  «  B(p)|u(x)  >  k) 

1  ;  {«  e  B(p)|u(x)  <  k) 

>«p 


(4.21)  /  IVulVdx  -  J  {f  (u  -  k)?2  -  f  u  C2  -  f.C*  2C(u  -  k)}dx 

\,P  *,P  1  1 

-  1  /.  wV( (u  -  k ) C 2 ) udx  -  /  (u  -  k)Vu?C2dx  . 

N.p  S.p 

The  Inequality  ab  <  «_1a2  ♦  «b2  and  the  facts  that  0  <  C  <  1  and  dlv  v  -0  Imply  that 
the  right  hand  side  of  (4.21)  Is  bounded  by 


2  l,0  *  (1  +  »)  I  2«  /*  l?u|Vdx 

AT  „  1-1  AT 

•t.P  ^t,P 

♦  jL  (u  -  k)2{1  ♦  (2  ♦  f)|Vd2}dx  ♦  X  !  IvlIVcKu  -  k)2dx 
*k.P  *k,P 

Choosing  J  •  1/2  and  observing  that 


!  *2dx  <  <pl  . 

Vp  2 

from  (4.21)  we  deduce  the  Inequality 


V 


(4.22) 


lV(u 


"‘■L.......  ♦  I, 

k,p 


|Vu|2CJd*  <  Yj(op)  2.« 


♦  If 


Kp' 


2 

2+k 


Yj(op) 


-1 


k.P 


|v|(u 


k)2dx 


where  tj  la  a  constant  depending  only  upon  the  data  and  la  Independent  of  v,  and  (op) 
la  aasuaed  to  be  leaa  than  or  equal  to  one. 

Proposition  3 

•  2  2 
There  exists  a  positive  number  s  depending  upon  the  data  and  the  [L  ((])]  norm  of 

♦ 

v  such  that  aithar 


aas  osc  a 


<  2 


•**/ 7+* 


B(R/2 ) 


or 

ess  osc  u  <  (1  *-2  ess  osc  u  . 

B( R/2 )  B(2R) 

This  proposition  is  proved  vie  a  aeries  of  laaeuts  which  make  use  of  the  inequality 

(4.22),  Setting 


U  ■  ess  sup  u,  v 
B(2R) 


ess  inf  u,  u  •  ess  osc  u  =  u  -  v 
B(2R)  B(2r) 


we  have,  obviously,  either 

(4.23)  I  (x  e  B(R)  t  u(x)  >  ||  -  |}  <  j  |  B(R>  | 

or 


(4.24)  |{x  €  B(R)  | u( x)  <  |l  -  |}|  <  -j  | B(R)  |  . 

It  is  assuMd  that  (4.23)  holds  with  the  arguments  being  similar  if  (4.24)  holds. 
Betting,  s  «  «, 


A 

s 


(4.23)  implies  that 


<4.25) 


I A  |  <  ~  WR2  »•  e  ■  . 


The  following  lean*  due  to  De  Giorgi  [5]  will  be  used. 


There  exists  a  universal  constant  8  such  that  for  every  pair  of  numbers  t,  k  such 
that  t  >  k, 

<i  -  k) I (x  e  B(R) |u(x)  >  t)| 1/2 


|{x  e  B(R)|u(x)  <  k}|  B,R)n[u>k)\[u>t) 


For  every  6Q  e  (0,1)  there  exists  an  integer  sQ  “  s0(v,60)  >  0  such  that  either 

w  <  2‘V/2~ 


I A  |  <  8n»R*  . 
s'  0 
0 


Applying  Iona  3  with  1  -  u  -  <#2  *  1  and  k  *  u  ”  »2  *,  s>1  yields. 


<4.26)  (^)  ,Vi|V2<_i5_r_j  |Vu|dx> 

•  V  8+1 

2 

By  virtue  of  (4.25)  |B(R)\Aal  >  vr  /2,  hence  an  application  of  HBlder's  inequality  to  the 
right  hand  side  of  (4.20)  results  in 


(4.27) 


Now  the  inequality 


'\+i' « 'Wi'!..  |7u' 


V\+i 


/  |Vu|2dx  <  /  |Vu| 2dx  »  I7(u  -  v  ♦  S-)4!2 

A  \ A  .  A  2* 


holds  and  this  last  tena  can  be  estimated  by  making  use  of  (4.22)  written  for 


k  •  li  -  #2  ,  p  •  2R  and  o  »  1/2  which  leads  to 


"<•  -  “  *  '-.>''1.,.,  <  vx*  - « * 


♦  v»*.  ' w/1~ 

*  U-=-,2R  B(  2R) 


Y,R_1  /  |v|(o  -  W  ♦  — — )+2d* 

*  atSek  «a® 


<  4»(~)2  +  Y  (4rS) 

2* 


2_ ,  1-2/2  ♦«  -1/»  -i 2 


♦  V1^)2  /  Ivldx 
2“  B(2R) 


since  (u  -  ii  ♦  «2  *)+  <  «2  *  on  B(2R).  By  HSlder's  Inequality 

R-1  /  |v|dx  <  R-’lvl  (v4R2)1/2  , 
B(2R)  2,Q 

so  that 


-  * ♦  ^+,2.B<R) 4  v’ *  ♦  u 


2k 

*-)2  ♦  y.r2+k 


which  leads  to,  upon  substitution  into  (4.27)  and  (4.28),  the  inequality 


(4.29) 


(prl'lv.'  <  (?>’lv> .  ♦  v'Tvv.' 


2k 

fin  ■> 2  .  „  _2+Ki 


If  RK/2+K  <  m2  0  and  Sg  >  s  than  (4.29)  iaplies  that 


(4.30) 


l*B+l  I  <  Y5(1  ♦  vm|\»|t1l,  s  -  1,2,  •  •  •  ,Sg  -  1 


adding  these  inequalities,  recalling  that  I  A, I  >  I A  I  and  |Aal  <  »R  /2, 

0 


obtain 


(S  -  1 )  |  A  |  <  Y5(1  ♦  V)wR  /2  . 

*0 


Renee  given  9^  e  (0,1)  to  prove  the  lsn  let 


(4.31) 


where  { r]  denotes  the  largest  integer  contained  less  or  equal  to  r. 
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we  have  fro*  (4.33)  that 


J+  l’u|2dx  -  /  |7(U  -  k  )*|2dx 

\  R  B(R  > 

n'Vl  n+1 


,4-”’  ‘  V4""  *  v,ll\../-2  *  K../  2"'«2  *  '< 

H  **  n  n  n  n 

wh*r®  Y6  *a  Independent  of  v  and  n.  Refining 


5n  “  K  .R  '/r2 
n  n 


(4.3S)  becomes 


„  _2  ( 1  -2/2+ic ) 

6n+t  <  V  (1  +  V>(^  +  5 - 5 -  «“ 


Provided  9.  <  1/w  we  have  that  0  <  S  <1  and  if 


♦  C  /2>  • 


then 

(4.36) 

where 


H  >  RK/2+K/2 


Sn+1  «  V^1  +  V)6T 


t  -  mini  4#  1  -  ~ — )  . 

On  the  basis  of  lenna  4.7  of  [10,  pe  66]  it  follows  that  provided 

«„  <  lT?(1  ♦  V))~1/T(24)'1/t2  , 
the  sequence  6n  converges  to  zero.  Hence  if  we  take 

<«•”>  90  -  I  (t7(1  ♦  V))-’/T(24)-,/t2} 

then  provided 

r2so  “  ,,  I  <  V*2  > 


(I) 
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either 


H  <  i  rk/2+k 


|A  iu  1  I  “  0  • 

V  -  — -  +  ~  H.R/2 

2  « 


Let  8n  be  chooen  as  in  (4.37).  By  lesma  4  there  exists  an  integer  sqOq'*)  such 


Proof  of  the  Proposition  3 

9o 


that  either 

(4.38) 
or 

(4.39) 


_  0_k/2+k 

(i)  <  2  R 


|A  |  <  8  ir‘ 

•o 


holds.  Tt>  prove  the  proposition  we  need  only  consider  the  second  possibility  and  show  that 
it  implies 
(4.40) 


-(sQ*1> 

ess  osc  u  <  (l  -  2  )  ess  osc  u 


B(R/2) 


B(2R) 


Prom  lemma  5,  (4.39)  implies  that  either 
(4.41) 


H  <  I  RK/2+K 


(4.42) 


)  A  .  |  -  0  . 

y  — ~  ♦  )«iV2 
2  0 


Suppose  (4.41)  holds.  By  the  definition  of  B, 


,  iii  1  k/2+*  .  iii  u  u 

ess  sup  u  <  y  -  -  +  —  R  <  y  -  —  +  - —  «  u  -  - - 

“o  2  *0  V1  V1 

2  2  2  w  2  u 


B(R) 

since  (4.38)  does  not  hold.  Hence 


ess  osc  u  <  ess  osc  u  <  «i  - 


B(R/2) 


B(R) 


.V1 


and  we  have  that  (4.40)  holds.  Otherwise  if  (4.42)  holds  then 


-23- 


i 


OJ  1 

ass  sup  u  <  U  -  -  ♦  T  H 

B( R/2 )  280 

a)  1  g  .  u  u) 

‘ 11  •  7T 5 

and  again  (4.40)  holds.  Consequently  the  number  s  claimed  by  the  proposition  Is  the 
Sg  given  by  (4.31)  and  (4.37). 

Lemma  6 

u  Is  Holder  continuous  in  S  with  exponent 

(4.43)  a  -  mln{^  *  —  >  In  1/o}  where  0-1-2  . 

Proof. 

For  every  t.  e  S  and  for  all  R  <  1  such  that  B(2R)  c  0,  set  u(R)  -  ass  osc  u. 

;  x/2tx  »<*> 

Me  have  shown  that  for  all  such  R,  either  ti»(R/4)  <  2  R  or  u(R/4)  <  ou(R) .  Then 

the  lean  is  an  issaediate  consequence  of  Isasa  4,8  of  t10,  p.  66]  and  in  particular 

(4.44)  osc  u  <  0  <  p  <  R 

B(P  ) 

where 

(4.45)  G  -  4amax{2Iul  _l  2*r’C^2+IC}  , 

Remark 

If  the  boundary  datum  ♦  in  Holder  continuous  with  *<skn  exponmnS  8  then  the 
demonstrated  H51der  continuity  of  u  carries  over  to  *  with  a  new  HBlder  exponent 
depending  on  V  and  B.  The  proof  of  this  follows  from  the  arguments  of  (10,  p.  90-95] 


with  the  modifications  indicated  above 


1  2 

Let  (lownti  of  the  quotient  space  W  '  (£J)/«  be  denoted  by  (n) .  The  quotient  non 

I  (nil  -  inf  mi 

netn]  w  '  (0) 

is  equivalent  to  ■  I  Vt»  1 1 2  q  for  *  tnl  -  Thus  for  any  u  6  L  (0), 

/  k,(u)VCVndx  [5), [nl  e  w1,2(n)/» 
a  c 

1  2 

defines  an  inner  product  for  which  W  *  (fl)/R  is  a  Hilbert  space, and 

/  k  (u)q(u,*)7nd*  v[n)  e  w1,2(a)/» 

Q  C 

defines  a  continuous  linear  functional  on  w1,2(a)/*.  Hence  the  Lax-Mil gram  theorem 
implies  that  there  exists  a  unique  (pi  e  H1' 2 ((!)/■  such  that 

(4.46)  p  e  Ip] ,  /  k  (u){7p  +  g(u,x)}Vndx  »  o  vn  e  [i]  e  w1,2(n)/R  . 

a 

The  pressure  p  =  p(u,e)  is  detenined  up  to  an  arbitrary  constant  by  (4.46)  whereas  the 
gradient  of  p,  and  hence  v(u,C)  -  -k£ (u) (7p{ u.r )  +  g(u,x)},  is  uniquely  determined.  It 
remains  only  to  obtain  the  Inequalities  (4.5)-(4.B).  Taking  n  ”p  in  (4.4)  yields  (4.5) 
and  hence  the  definition  of  v(u,c)  gives  (4.6).  Inequality  (4.7)  is  proved  by  writing 
(4.4)  for  u^  and  Uj,  differencing  and  choosing  n  *  p(u^,e)  -  p(u2>e).  Finally  (4.8) 
is  a  consequence  of  (4.7)  and  the  definition  of  v(u,e). 


5.  PROOF  OP  PROPOSITIONS  1  AMD  2 


•  )  2 

Fix  e  >  0  and  u  e  L  (ft).  Let  p(u,c)  e  W  '  (0)  solve 

(5.1)  -  div  k£  (u)  (Vp(u,e)  ♦  g(u,x)]  -  0  in  V’ {SI) 

-k£(u)  (Vp(u,e)  +  g(u,x)]N  -  0 

♦ 

and  v(u,e)  be  uniquely  defined  by 

(5.2)  v(u,e)  -  -kE(u) (Vp(u,e)  ♦  g(u,x)}  . 

Set  F(u)  e  W1,2(8)  to  be  the  unique  solution  of 

(5.3)  -AP(u)  +  Xv(u,e)Vp(u)  -  fQ  ♦  t  ^  in  £>•  (0) 

P(u)  “  ♦  on  3ft  . 

Obviously  (5.1)-(5.3)  are  meant  in  the  weak  sense  made  precise  earlier. 

Lemma  7 

(i)  F  -  l”(«)  ♦  L** {SI)  is  well  defined. 

(ii)  There  exists  a  constant  M  depending  only  upon  the  data  in  [Aq]  - (A^)  such 
that  if  B(M)  is  the  ball  of  radius  H  in  L  (ft),  then 

F  s  B(M)  ♦  B(M)  . 

(iii)  F  is  completely  continuous. 

Proof.  Lemmas  1  and  2  imply  that 

.F,u).. ^  <  T , <  1  ♦  2k0{C|fl|1/2.U.^n  ♦  .,0.2#a})  H  7,  ♦  Y^lul*  fl  . 

Let  M  be  the  solution  of 


Y,  ♦  -  N  >  0  . 


Such  an  M  exists  since  t  6(0,1),  and  it  can  be  determined  a  priori  in  terms  of  the  data 

only.  Therefore  if  u  e  B(M)  we  have  F(u)  6  B(M).  This  proves  (ii).  Statement  (i)  is  a 

♦ 

consequence  of  the  uniqueness  of  v(u,c)  and  F(u)  defined  by  (5.2)  and  (5.3). 

By  lemma  1  F(u)  is  Holder  continuous  on  compact  subsets  of  ft  with  constants 


depending  upon  Ivl^  ft  *  u  6  B(M)  then  by  (4.6) 


"  ■  t 

,v,2,n  *  T1  *  V 


where  y  and  y2  depend  only  upon  the  data.  Therefore  the  HSlder  constants  of  Flu) 


can  be  a  priori  determined  depending  only  on  the  data.  Consequently  F  is  compact.  In 


order  to  show  continuity  of  F,  let  (u^}  be  a  sequence  in  B(M)  such  that 
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Uj  ♦  uQ  8  BIN)  uniformly  In  !)•  Writing  (S.l)  in  the  Mik  fora  for  Uj  and  uQ  we 
obtain,  by  differencing,  the  equation 

/  <V(r(Uj)  -  r( uQ ) )  -  l»(Uj,t)(r(Uj)  -  P(u0))}Vndx 

-  -  /  {w(u  .£)  -  v(u  ,t)}vr(u  )n<bt  vn  e  w1,2<fl>  . 
o’  0  0 


We  take  n  ■  F(  u  ^ )  -  f(ug)  and  observe  that  the  second  Integral  on  the  left  hand  aide 
vanishes  for  this  choice  of  n,  since  div  v  ■  0  in  D'(Q).  Thus  we  obtain 

/  IVWu^J  -  r<u0Jf2d*  <  GOi.ivriu^lj^JivtUj.el  -  »<«0«e>'2f8 
In  view  of  (4.7)  and  (4.B),  the  following  inequality  holds 


iV(r(Uj)  -  r(u0))i2#fl  «  G(ii.ivr(u0)i2#Q,d.ta.e).(iuj  -  u0i_  fl  ♦  i|oj  -  c0H.>fl} 

where  -  gtu^x).  Letting  j  ♦  «  wo  see  that,  since  g(»,*>  is  continuous  by  [hjl , 

f(Uj)  ♦  MUg)  strongly  in  tf*'2^) 

However  the  uniform  boundedness  and  HSlder  continuity  of  (FCu^))  implies  for  every 

subsequence  (f(u',)}  the  existence  of  a  subsequence  u  ,  for  which 
3  3k 

F(u., )  ♦  F  in  l" (0)  . 

3k 

Hence  F  «  F( uQ )  and  for  the  entire  sequence 

F(uj)  ♦  F(u0)  in  L  (0)  . 

F  is  continuous  and  the  lemma  is  proved. 

•ft 

By  the  Leray-Schauder  fixed  point  theorem  f  has  a  fixed  point  in  L  (0).  Therefore 
we  can  conclude  that  Vt  >  0  problem  (2.2)-(2. 5)  has  a  solution  and  Propositions  1  and  2 
then  follow. 


-27 


gjgjrt 


In 


aull> 


tht  cbm  t  m  1#  there  exists  en  M  >  0  solving  (5.4)  when 

IkgCini1^2  <  1.  Thus  we  hsve  existence  for  t  -  1  when  the  dsts  is  sufficiently 
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